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We investigate the ground state of the Kondo necklace model on geometrically-frustrated lattices by 
the variational Monte Carlo simulation. To explore the possibility of a partially-ordered phase, we employ 
an extension of the Yosida-type wave function as a variational state, which can describe a coexistence of 
spin-singlet formation due to the Kondo coupling and magnetic ordering by the Ruderman-Kittel-Kasuya- 
Yosida interaction. We show the benchmark of the numerical simulation to demonstrate the high precision 
brought by the optimization of a large number of variational parameters. We discuss the ground-state phase 
diagram for the model on the kagome lattice in comparison with that for the triangular-lattice case. 
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1. Introduction 

Competition between different electronic phases is a key to 
understand the variety of fascinating phenomena in strongly- 
correlated electron systems. An interesting situation shows up 
when different phases compete with each other at zero tem- 
perature (T) to result in a quantum critical point (QCP). QCP 
is a source of many exotic properties, such as a non-Fermi- 
liquid behavior and some emergent phases. 

A typical example of QCP and related phenomena is found 
in the Kondo lattice systems such as /-electron systems.'' In 
the Kondo lattice systems, competition originates in two dif- 
ferent types of interactions; one is the Kondo coupling be- 
tween conduction electrons and localized moments, and the 
other is the Ruderman-Kittel-Kasuya-Yosida (RKKY) inter- 
action which is an effective magnetic coupling between local- 
ized moments mediated by conduction electrons. The former 
induces the formation of the Kondo spin-singlet,^* whereas 
the latter tends to stabilize a magnetic long-range ordering.^' 
The competition leads to QCP between a Fermi liquid state 
and a magnetically-ordered (MO) state, which has been ar- 
gued as a key concept for understanding rich behaviors in the 
Kondo lattice systems.''* 

Geometrical frustration is another important concept in 
strongly-correlated electron systems. Frustration suppresses 
conventional magnetic orderings and opens up the possibil- 
ity to have exotic phases, such as fluctuation-induced phases 
by the order-from-disorder mechanism.^' One of such intrigu- 
ing phenomena is a partial order:^' The system is separated 
into magnetic and nonmagnetic sublattices to relieve the frus- 
tration; namely, a magnetic ordering appears on a sublattice 
of the system with leaving the remaining sublattice nonmag- 
netic. The partial ordered states are in most cases stabilized 
by thermal fluctuations, i.e., by a gain of the entropy. 

Recently, a different type of partial order was explored in 
the Kondo lattice system by the authors.'* The partial order 
is not entropic-driven but a new quantum state stabilized at 
zero T under geometrical frustration in the competing re- 
gion between the Kondo coupling and the RKKY interaction. 
The partial ordered phase consists of a magnetically-ordered 
unfrustrated network in a sublattice and nonmagnetic spin- 



singlet sites with Kondo screening in the remaining sublattice. 
We call this state the partial Kondo screening (PKS) state. The 
phase diagram was studied by the variational Monte Carlo 
(VMC) method, and it was shown that PKS is further stabi- 
lized by quantum fluctuations and the spin anisotropy. 

In this contribution, we report the numerical details of the 
VMC calculations. By showing the extended benchmark re- 
sults with detailed conditions of the calculations, we demon- 
strate that VMC gives the ground-state energy in high pre- 
cision, typically within 1-3% relative errors. We also show 
a possible ground-state phase diagram for the Kondo neck- 
lace model (KNM) on the two-dimensional kagome lattice. 
We discuss the stability of PKS in comparison with the re- 
sults for the triangular-lattice case. 

2. Model and Method 

In the present study, we consider KNM which has been 
studied as a mimic of the half-filling case of the Kondo lat- 
tice model."'* The Hamiltonian is given by 
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where the first term describes the AF interaction between con- 
duction electron spins r^ and the second term represents the 
Kondo coupling between Ti and the localized spin Si. The for- 
mer tends to stabilize a magnetic ordering, and the latter fa- 
vors the Kondo spin-liquid (KSL) state in which all the sites 
participate in the local Kondo singlet formation. Hence, in the 
unfrustrated cases, this model describes QCP between MO 
and KSL, instead of MO and a Fermi liquid state, because of 
the half filling."'* For simplicity, we consider S = 1/2 spins 
for both fi and Si. In Eq. (1), we extend the model by adding 
the last term, the AF Ising interaction between the localized 
spins {S^ is the z component of Si), in order to incorporate 
the magnetic anisotropy often seen in real materials. We con- 
sider the model on two different frustrated lattice structures 
in two dimensions, the triangular lattice and the kagome lat- 
tice, shown in Fig. 1. The sums (ij) in Eq. (1) are taken over 
the nearest-neighbor pairs on these lattices. Hereafter we take 
W^ = 1 as the energy unit. 

We study the ground state of the model in Eq. (1) by the 
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Fig. 1. (Color online) Schematic pictures of three-sublattice orderings on 
(a) the triangular lattice and (b), (c) the kagome lattice, (b) is the q = 
ordering and (c) is the \/3 X \/3 ordering. 12- and 24-site clusters used in 
the present calculations are also shown.'" 



VMC method. To explore the PKS state in which nonmag- 
netic spin singlet and magnetic ordering coexist, we here con- 
sider an extension of the Yosida-type wave function^' as the 
variational wave function. The functional form is given by 

IV> = PG/:^="/:^="l0pah-). (2) 

Here |(/)pair) is a generalized BCS wave function given by 
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where /^m are the variational parameters, N is the number 
of lattice sites, and |0) is a vacuum. The pair creations with 
fermion operators a\^ include both conduction and localized 
spin degrees of freedom (the sums run over the lattice sites 
with distinguishing them); the electron numbers are fixed to 
be one in each degrees of freedom by the Gutzwiller projec- 
tion operator Pq in Eq. (2). £^=" and C^^^ are the quantum- 
number projection operators for the total spin singlet and the 
total momentum zero, respectively. We optimize a large num- 
ber of variational parameters by using the stochastic recon- 
figuration''^ and enforce the quantum-number projections, by 
following the method in Ref. 10. 

3. Benchmark 

In this section, we show the benchmark results of our nu- 
merical calculations. In the following, we explore the ground 
states with three-sublattice ordering in two different ways. 
One is to impose the momentum projection C^ ^^ only for the 
same sublattices with taking account of all possible /;„. The 
other is to impose the three-sublattice translational symmetry 
on Jim, which largely reduces the number of Jim required in 
the calculations. We take the 12- and 24-site clusters shown 
in Fig. 1 under the periodic boundary conditions compatible 
with the three-sublattice orderings. We apply the spin-singlet 
projection C^^^ only for the case with 1^=0 hereafter. We 
typically perform 300-1000 stochastic reconfiguration steps 
with 1600-6000 MC samplings in the following calculations. 
The benchmark is shown for the parameter range where the 
three-sublattice orders (the MO and PKS states) are obtained. 

First we present the results for the triangular-lattice case. 
Table I shows the total ground-state energy for the 12-site 
system obtained by VMC in comparison with the result by 
the exact diagonalization. The VMC results show the ground 
state energy by the first method, i.e., by the full optimization 
of fim with £^^". As shown in the table, VMC provides the 
energy very close to the exact value in a wide range of param- 



Table I. Ground-state energy for KNM on the triangular lattice with A^ = 
12. Comparison between the VMC results and the exact diagonalization 
results are shown with relative en'ors. The numbers in the parentheses indi- 
cates the statistical errors in the last digit. The upper table is for the results 
aX Iz = and the lower at Iz = 0.4. 



J 



0.2 



0.4 



0.6 



0.8 



VMC 
exact 
eiTor 

J 


-7.38(1) 
-7.446 
0.009(1) 

0.2 


-7.98(1) 
-8.061 
0.010(1) 

0.4 


-8.93(1) 
-9.025 
0.011(1) 

0.6 


-10.16(1) 
-10.232 
0.007(1) 

0.8 


VMC 
exact 
error 


-8.49(1) 
-8.602 
0.013(1) 


-8.92(1) 
-9.103 
0.020(1) 


-9.72(1) 
-9.957 
0.024(1) 


-10.77(1) 
-11.059 
0.026(1) 



Table II. Ground-state energy for KNM on the kagome lattice with N 
12. The upper table is for Iz = 0, while the lower for Iz = 0.4. 

J 0.14 0.16 0.18 0.20 



VMC 
exact 
error 

J 



-5.55(1) -5.58(1) 
-5.578 -5.642 

0.005(2) 0.011(2) 



-5.67(1) 


-5.76(1) 


-5.717 


-5.800 


0.008(2) 


0.007(2) 



VMC 
exact 
error 



0.14 



0.16 



0.18 



0.20 



-6.30(1) -6.33(1) 
-6.349 -6.403 

0.008(2) 0.011(2) 



-6.36(1) -6.40(1) 
-6.469 -6.544 

0.017(2) 0.022(2) 



eters J and Iz- The relative errors are typically ^ 1-3%. 

Table II presents similar comparisons for the 12-site 
kagome-lattice models. Note that the system with A^ = 12 un- 
der the periodic boundary conditions accommodates only the 
q — type three-sublattice ordering, as shown in Fig. 1(b). 
As presented in the table, VMC gives the ground-state energy 
in high accuracy also in this kagome-lattice case, typically 
within 1-2% relative errors. 

Finally, we compare the two methods, namely, the full opti- 
mization of fim with £^=° and the optimization of a reduced 
set of fim with three-sublattice periodicity. The results for the 
24-site triangular-lattice models are shown in Table III. For 
the 24 site, the number of the variational parameters /^ is 
2304 for the full optimization, whereas it is reduced to 288 
when the three-sublattice periodicity is imposed. The results 
show that the reduced set still gives good precision, '^^ in par- 
ticular, for larger J where spatial correlations become local. 
This encourages us to extend the calculations for larger clus- 
ter sizes by using the reduced set of /;„ which greatly reduces 
the CPU cost compared to the full optimization. The extension 
will be reported elsewhere. 

4. Results for the Kagome Lattice Case 

Here we present the VMC results for KNM on the kagome 
lattice, in comparison with those for the triangular-lattice case 
previously reported in Ref. 7. In the kagome case, even when 
considering only the three-sublattice ordering, there are dif- 
ferent ordering patterns, such as (f = or a/3 x \/3 [Figs. 1 (b) 
and (c)]. This is because of the comer-sharing network of tri- 
angles. Since it is difficult to analyze the ground state with 
allowing all such possibilities in VMC with increasing the 
system size N, here we limit ourselves to a small size clus- 
ter iV — 12 under the periodic boundary conditions, which 
accommodates only the solutions compatible with the q = 
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Table 111. VMC ground-state energy for KNM on the triangular lattice with 
Af = 24 sites. The results show the comparison between the optimization 
of all /jm and of only three-sublattice symmetry-allowed fim- Smaller 
errors in the former cases are owing to the quantum projection C^~^. 
The upper table shows the results at I^ = 0, while the lower at I^ = 0.4. 



J 


0.2 


0.4 


0.6 


0.8 


all fhn 


-14.20(7) 


-15.15(4) 


-17.03(4) 


-20.04(8) 


3-sub /i„ 


-12.2(3) 


-15.0(2) 


-17.1(3) 


-19.1(3) 


J 


0.2 


0.4 


0.6 


0.8 



all flm 
3-sub /i„j 



-16.62(3) -17.60(4) 
-14.6(3) -16.8(1) 



-19.31(3) -21.49(7) 
-18.6(5) -20.9(2) 




Fig. 2. (Color online) Phase diagram deduced from the VMC results for 
KNM on the kagome lattice with Af = 12 sites. The data for the triangular- 
lattice case with N = 12 are taken from Ref. 7. Schematic pictures for 
the spin state in each phase are shown, in which the grey (black) an'ows 
represent r; (Si) and the circles denote the Kondo singlets. 



type order. The following calculations are done by the full op- 
timization of all /;,„ with imposing £^^ ^'^ for the same sub- 
lattices. We map out the ground-state phase diagram from the 
behaviors of the onsite spin correlations (r^ -Si) as well as the 
intersite ones {Si ■ Sj), following the manner in our previous 
study for the triangular-lattice models.^* 

Figure 2 shows the results. For comparison, the results for 
KNM on the triangular lattice with TV = 12 are also pre- 
sented. Note that the horizontal axis denotes J renormal- 
ized by the coordination number z (the number of nearest 
neighbors) to facilitate the comparison. We obtain a simi- 
lar sequence of three different phases as in the triangular- 
lattice case, i.e., MO, PKS, and KSL phases, with increasing 
J. However, we observe that the PKS phase is largely sup- 
pressed in the kagome-lattice case, and that the KSL phase 
is relatively stabilized in the large J region. In particular, 
we do not see a clear evidence of PKS at I^ — 0, whereas 
the triangular-lattice case exhibits a finite window of J for 
PKS down to Iz = 0. The results suggest that PKS is rela- 
tively unstable in the kagome-lattice model compared to the 
triangular-lattice model. The weak tendency to PKS is rea- 
sonable since the magnetic sublattice in PKS has a relatively 
weak network in the kagome-lattice case compared to the 
triangular-lattice case: it forms a two-dimensional honeycomb 
network in the triangular-lattice case,^' whereas it consists of 
one-dimensional chains in the q = ordering in the kagome- 
lattice case. The tendency is also consistent with the observa- 



tion in the limit of /^ ^ J (^ W) where the system can be 
mapped onto an effective transverse-field Ising model:^' For 
the transverse-field Ising model, a three-sublattice partial or- 
der is suggested to be stabilized in the triangular-lattice case, 
whereas it is not in the kagome case.'-'* 

Even though we believe that the results capture the correct 
tendency, it is unclear whether the PKS state survives in the 
kagome-lattice model in the thermodynamic limit iV — > cx). 
Especially, it will be important to take account of all other 
possibilities of the three-sublattice ordering patterns. The 
most exotic possibility is a spin-liquid state composed of the 
PKS three-site units (two sites bear magnetic moments, while 
the remaining one participates in the Kondo-singlet forma- 
tion): In the kagome lattice, the triangle units are connected 
with sharing only their corners, and hence, it is possible to 
have a magnetically-disordered state by connecting the PKS 
three-site units. In the extreme case, such state bears macro- 
scopic degeneracy. It is interesting to ask whether the degener- 
acy is lifted through quantum order-from-disorder. These in- 
teresting issues are left for future study. 

5. Summary 

We have investigated the ground state of frustrated Kondo 
necklace model by the variational Monte Carlo calculations. 
We have shown the detailed benchmark for both the triangular 
and kagome lattice cases, which demonstrates that our sim- 
ulation employing a generalized Yosida-type wave function 
gives highly-precise ground-state energy in a wide parame- 
ter range. We have also deduced the phase diagram for the 
kagome-lattice model. The comparison with the triangular- 
lattice case revealed the reduced tendency toward the partial 
Kondo screening in the kagome-lattice case. 
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